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'— ' Abstract. 

, Wc present a simple kinematic model of a non-equilibrium steady state device, 

^ which can operate either as a heat engine or as a refrigerator. The model is composed 

^"^ of two or more scattering channels where the motion is fully described by deterministic 

^ 1^ classical dynamics, which connect a pair of stochastic (infinite) heat and particle baths 

1/-^ at unequal temperatures. We discuss precise kinematic conditions under which our 

I I model may approach Carnot's optimal efficiency in different situations. 
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1. Introduction 

It has been recognized in recent years tliat an efficient energy conversion and 
transformation between heat and work on all scales, in particular at the nanoscale, 
poses one of the main challenges for future social needs. One of the most promising 
technologies is a thermoelectric circuit which can simply convert a heat flux from a hot 
to cold place into an electric power, or conversely, using the electric power to pump 
heat from cold to hot place pj. However, although such devices have been known for 
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a long time, and even their technological applications have been diverse for several 
decades, their thermodynamic efficiency remains quite low compared to compressor 
based refrigerators or fuel burning internal combustion engines. One of the main reasons 
for such low efficiency is a rather poor theoretical understanding for the relatively low 
figure of merit, the so called ZT in all known materials at room temperature. 

Recently, a dynamical system approach has been proposed [6|, [11] in an attempt 
to better understand the mechanism which could lead to an increase of ZT and 
thus of thermoelectic efficiency. Along these lines we propose here a mathematically 
simple, general model of a steady state thermoelectric device, which can be analyzed 
theoretically within the framework of classical statistical mechanics. The model is 
composed of two or more classical two-port scattering channels where the motion (of non- 
interacting charged particles) is fully deterministic and energy conserving, and which 
connect two stochastic, infinite heat and particle baths at different temperatures. The 
motion in the longitudinal direction connecting the two baths is quasi-free except for 
point scatterers in the middle of the channels, which have certain simple properties. 
Namely the transmission only depends on the (kinetic) energy of the particle. It is 
argued that the dependence on energy of the transmission (scattering) function is crucial 
in order to obtain non-vanishing efficiency. The particles are only subjected to a bias 
voltage that is used to input or output the work from the system. 

The model we discuss here, even quite abstract, is the most general case of a self- 
contained (compact) classical heat engine model, composed of two thermo-chemical 
baths with non-interacting particles. Its simplicity allows theoretical analysis and 
therefore a better understanding of how a thermo-electric engine operates and the 
mechanism underlying the behavior of ZT. Indeed in our paper we develop an easy- 
to-follow formalism in order to discuss particle and heat currents inside the model as 
well as its power and thermoelectric efficiency. The model's simplicity indicates that 
engines based on its principles could be practically built and so the model may be 
important for future design of heat engines. 

More precisely we derive analytical expressions for the engine's efficiency, power and 
the parameters that optimize efficiency (section 2). In the linear temperature difference 
regime (section 3) in particular, there are cases like window-type scattering functions 
(section 5) that can be worked out exactly. It is particularly interesting that in the linear 
temperature difference regime, refrigeration and power generation modes of operation 
can be treated symmetrically. Using our formalism, it is possible to discuss the precise 
conditions under which the model attains Carnot's efficiency. We show in section 4, that 
a generalized figure of merit is essentially given as a ratio of the variance of Seebeck 
coefficients for various channels divided by the average variance of the transmitted 
particle's energy. It is found that narrow energy filtering is a good mechanism for 
increasing thermoelectric efficiency of the engines based on our classical model. Similar 
effect was previously observed in quantum models, see e.g. p]. In section 6 we worked 
out semi-analytically and numerically a slightly different model of thermoelectric engine 
where the scattering is provided by a homogeneous transverse magnetic field. In this 
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setup, power and efficiency turn out to be tunable via the external magnetic field; 
however the maximal efficiency remains rather low. 

Some more technical aspects of our formalism which are not essential for the 
main results are presented in two appendices. In appendix A we present consistent 
construction of classical scattering matrices and their concatenation for structured 
scattering channels. While most of the work in this paper refers to the steady-state 
operation of heat-engine models, we dedicate appendix B to the discussion of relaxation 
times to the steady state. In particular we show that non-vanishing bias voltages in all 
channels render strictly finite relaxation times. 

As discussed in the following sections, our formalism allows an efficient numerical as 
well as analytical treatment of concrete not trivial cases of classical-mechanical models 
of heat engines and helps us to identify the main transport properties controlling power 
generation and efficiency. 

2. Construction of the heat engine or refrigerator 

The heat engine or refrigerator is constructed of two d- dimensional thermo-chemical 
baths connected by two or more channels as depicted in figure [1} The particles in the 
system have d degrees of freedom, equal positive charge e and equal mass m. We use 
units in which the charge e, mass m and Boltzmann constant k^ are equal to unity, 
m = e = /cb = 1- Along the direction of the channels we introduce bias electrical 
fields Ei, which are described by bias voltages f/j between the left and right baths. 
These voltages are used to extract work (or power) from the system in the case of 
power producing heat engine or to supply the work in the case of the refrigerator. 
The Hamiltonian H oi a particle with momentum p = (pi, . . . ,pd) G M'^ and position 
q = {qi, . . . ,qd) G M'^ inside the ith channel is given by 

H{q,p) = W{p) + U,q^, W{p) = ^\p\'. 

The coordinate axis gi runs longitudinally along the channel with units chosen so that 

qi G [0, 1], where gi = and gi = 1 mark the left and the right end of the channel, 

respectively. Let there be n channels connecting the left bath side {u = L) with the right 

bath side [u = R) via the holes Fj C M'^"^ perpendicular (transverse) to the direction 

of the channel of "area" Ai = volFj. For convenience we assume that the holes through 

which the channel i is attached to the left and the right bath are represented by identical 

sets Fj (analysis could easily be generalized if necessary). The total area of the holes 

is A = Yl^=i ^i- The ratio Vi = Ai/A represents the probability that a particle effused 

from the bath goes into i-th channel, from either side. The dynamics in each channel is 

reversible and energy conserving. The transmission function for a particle from t/-side 

over ith channel is defined as 

particle is effused at x G Fj with momentum p ^M.'^ from side 
(x,p, Ui) = \ ^ ^ ^'^'^ transmitted through channel i with bias voltage Ui (1) 

: particle is refiected back into the same bath 
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Figure 1. Schematic figure of the heat engine. 



A scheme for explicit construction of the transmission functions is discussed in appendix 
A. In the left (right) bath the particles are at temperature Tl (Tr) and chemical potential 
/^L (/^r) and are effused into i-th. channels with the rate rj7L (rj7R). The effusion rates 
7l and 7r are connected with the chemical potential /Xy and reciprocal temperature 13^ 
via the formula 

/3/i^ = log {ci3^A^-i^^ , z/ G L, R , 

where C is a constant depending only on the mass of the particles. The particles effused 
from i/-bath into channels have the momentum j9 G M'^ distributed according to the 
Boltzmann distribution 



Pip.iyJ) 



' V27r 



bilexp ( --/3|pp ) 9{api) , 



(2) 



with unit step function 0{x) = (0 : x < 0; 1 : otherwise). Let us now introduce the 
probability t^^i that a particle effused from side v is transmitted over i-th channel, and 
the average particle kinetic energy gj,j, which are computed as 

{tu,i,qu,i) = X / dpP(p,cr,„/3j,)(l,iy(p)) / dxT^^i{x,p,Ui) , (3) 

A jRd JTi 

with (Tl = 1, cs^ = —1. Notice that at equal temperatures Tl = Tr and zero fields 
f/j = we have tL,i = tR,i and gL,i = ?r,«- The particle currents jp^ within i-th channel 
and the heat currents jq,i\u [Sj exchanged between the i/-side bath and i-th channel are 
given by 

jp,i = rii'jLtLA - 7R%,i) , (4) 

jq,i\-L = ri (7LgL,i - iRiqRA + tK,iUi)) , (5) 

iq,i\K = ri (7L(gL,i - tL,iUi) - 7RgR,i) . (6) 

In the stationary state, the engine model should obey two basic principles (compatible, 
or even equivalent to the 2nd law): 

• The net particle current from one bath to another is zero X]r=i^p.« ~ '-'• 

• If temperatures of the baths are equal, then the particle current in each channel is 



zero. 
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Assuming that the average injection rate 7 = ^(7l + 7r) is known (say it is an input 
parameter of the model) the two above principles yield the relation 

7r_-,_._7l_ a -Ar 

7 7 ^L + ^R ' 

which uniquely determine the injection rates 7l and 7r. After introducing auxiliary 
variables 



(^wAjy, Oj/, Uj/j / j^i\J^u,iy'^u,i^iyQu,i) ; 



1=1 



the stationary particle currents in the channels read 

^P'« ^ ~A , A~ [^L.i-^R ^ ^R,i-4,L] 

•^L + v^R 

and the power produced by the particle flow is given by 

n 

P = Y. U^h. = . T\ [^L^R - i^RA] . 
^ A.I, + >Ir 

1=1 

A non-vanishing power is produced, if and only ii T> = Bi^Ayi — B^Ai^ is non-zero. 
Similarly, we define the total stationary heat current Q^ = Y17=i3gAi' '^^ both sides of 
the channels 

Ql = , ^/, [Cl^r - CrA + BrAl] , (7) 

-^L + -Ar 

Qr = , ^\ [ClAk - CrA + Bi^An] ■ (8) 

vAl + vAr 

Notice the similar expression for the total heat current on the left and right side, which 
become identical in the linear temperature-difference regime. We can think of the 
apparatus as a power producing heat engine or as a refrigerator. In the heat engine 
regime, the heat is transferred from the hotter bath u to the colder bath in such a way 
that the particles in the channels travel against the bias voltage and generate a power 
P > 0. In the case of the refrigerator regime, the heat is pumped from the colder bath 
of the side z/ by appropriately setting the potentials in the channels and thereby a power 
P < is used. By introducing a binary index 

J +1 : in the case of the power producing heat engine mode 
I —1 : in the case of the refrigerator mode 

we can write the efficiency rj for both cases as 

""l^^j <*»o.= (^j— j . (9) 

with (TjyQi/ > 0. The efficiency rj of the heat engine and refrigerator is bounded from 
above by the efficiency of the ideal Carnot cycle r^camot- For a given setup of scatterers 
the efficiency can be optimized by choosing appropriate potentials. In general the 
optimization can be done only numerically, because the expressions for t,y^i and q^^i are 
usually nonlinear functions of the potentials. Only in the linear temperature-difference 
regime, discussed in the next section, the optimization can be done analytically. 
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3. Linear temperature-difference regime 

In the previous section we have discussed the general framework of our heat engine 
model. In the following we elaborate on the linear temperature-difference regime, 
where the temperature difference 6T = Tr — Tl is much smaller than the average 
temperature T = |(Tl + Tr). Consequently, the difference in injection rates ^7 = 
7r ~ 7l and in electrical potentials Ui on both sides are also small. These conditions 
permit a perturbative treatment of energy and particle currents and enable an analytic 
optimization of the efficiency. 

The idea is to expand the particle Q and the heat currents ([s]) (|6]), which depend 
on transport coefficients ts,i and qs^i ^ and injection rates 7^, in terms of ST, S'j and 
Ui- First we write the injection rates from both sides 

lu = l + ^(^Jl, z/e{L,R}. 



using the difference in injection rates ^7 reading 



57 = 7 



/i — ) 5/3 + /35/i 



(10) 



where we introduce the average and the difference of chemical potential /Z = (/iL +/^r)/2 
and 5/i = /iR — /^l, respectively, and the reciprocal average temperature (3 = T . 
Next we expand the transport coefficients tgi and g^.j linearly in /3 and f/j around the 
equilibrium point (3 = (3, Ui = 0: 

tu,i = tu,i\o + dj3t^^i\o 5(3 + ayduitv,i\o Ui (11) 

qu,i = qu,i\o + df3q„,i\o 5(3 + (JyduSvAo Ui ■ (12) 

Then, by using the transport coefficients at the equilibrium point 

{tf\qf\kf)=]^! dpP{p,lMhW{p),W{py) f dxT,,{x,p,0) , (13) 

we may write the temperature expansion coefficients as 

d+l jo) (0) ^ ^ I _ l±l^(o) r,(o) 

Ol3tu,i\0 - -5 t- - g ■ , C/3gj,,j|o - - Qi - Ki . 

The latter depend only on the mean reciprocal temperature (3. For the time being 
we assume that the channels are straight cylindrical leads with a local scatterer in the 
middle (i.e. located in the center of each channel) and that the scattering process 
depends only on the kinetic energy, i.e. 

T,,i{x,p,Ui) := ^i (w{p) - laMi] (Ipj -max{0,aM)] , (14) 



2 *^ 7 V2 

where ■?/' is a transmission function at fixed energy of the scatterer inside the channel. 
We refer to ip as the on-shell transmission function. It is possible to show that 
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and the transport coefficients at equilibrium point simplify to 

d+l 



(0) „(0) , (0) 



/3 



(¥)' 



dW W^e-^'^^i{W){l, W, W 



(16) 



The relations (15) are valid for scatterers of more general type even though it might 
be difficult to give a rigorous proof. In order to simplify the following discussion we 
introduce auxiliary constants 



.(0) 



.(0) 



U = ll^rit] ' , qi = -f(3riql , ki 

and the relative reciprocal temperature difference 



ifirik 



(0) 



(17) 



Taking into account the linear expansion of transport coefficients (15) and injection 
rates (10) we can write the particle currents jp^i (|4| and the heat currents jq^i ([s]) as 

3p,^= -ti{Ui + 6fi)+^qi, (18) 

Jq,^= -qi{Ui + 6fi)+^ki. (19) 

where we use the fact that the heat current on the left and right side are equal, 
jq,i '■= Jq,i|L = jqA'Ry ^^ ^^^ cousldered regime. Without loss of generality, the mean 
value of the chemical potential JI is set to zero, as only differences are important. The 
expression for the currents can be given in the matrix form 



] p,i 
Jg,i 



(20) 



qi h -{Ui + 6fi) 

showing that they indeed satisfy the Onsager reciprocity relations [3]. It is convenient 
to express qi and ki in terms of positive quantities Si = qi/ti and Ki = ki — tiSf and 
introduce n-component column vectors 



"i^i 



t = iti\u, S = [S,] 



i}i=l ' 



K=[K, 



«Ji=l ' 



f/=[f/.]iLi, i = [i,...,ir 



The stationary condition of zero net current ^"^^ jpA = determines the difference 6fj, 
of chemical potentials 



Sfi 



I't 



By using this result we write the stationary power P = J^iJpA^^i and the stationary 

as 



net heat current Q = ^ijq,i 

P={^S- UfMU , g = e [S^MS + (l^K) {ft 



MS . 



where the matrix M = (l^t) diag(t) — tt^. Notice that P and Q are invariant under the 
transformation U ^^ U + const. The efficiency ?7 (^ of the heat engine or refrigerator 
can be optimized w.r.t. to the field vector U by solving the set of equations drj/dU = 0. 
Let us introduce a quantity 

S^MS , , 

y= , ^ ,, ^ , >0 21 

l^fsT l^t ^ ' 
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called the merit of efficiency. Then the optimal efficiency rj* of the heat engine or the 
refrigerator is equal to 

yi + y-i .„„^ 

and the corresponding vector of the bias voltages U* and the output power P* are 
written as 

u* = e^^ (yrr^- c) (^ - ^(i^^)i) ^ (23) 

P* = ~ (j) {S^MS)^/lT^ (2v/rT^ - C(2 + y)) , (24) 

Q*=C-iS^MS),/l+^, (25) 

1/ 

where we use for Carnot efficiency r/camot = C^'' ^- The relative efficiency rj /rjcamot of 
the heat engine or refrigerator is a monotonically increasing function of y. In the limit 
y — > 00 we reach the Carnot cycle efficiency. In the limit of two channels ra = 2 we 
obtain the expression already derived in [SI [ID] • In this case the merit of efficiency reads 

(5i - S2Y 



4. Optimally efficient engines 

The efficiency of a steady state device that is working as a heat engine or refrigerator 
can be far from maximal, Carnot cycle efficiency. Here we describe a general strategy 
to increase efficiency which is based on our model in the linear temperature difference 
regime. 

Let us consider a ra-channel heat engine or refrigerator with scattering depending 
only on the kinetic energy T^^i{x,p, 0) = ipiiW^p)), at reciprocal average temperature (3. 
We want to determine which properties of the scatterers may lead to maximal efficiency. 
By taking into account the transmission function and the distribution of momenta by 
the effusion process ^ we know that the particle's energy W transmitted over ith 
channel has the probability density 

_d+i 

with average Wi = Si and variance Vi := W^ — Wi = Ki/ti. It is interesting to notice 
that in our notation the figure of merit ZT [6] of the ith channel is given simply by 

In order to further simplify the discussion we assume that the bath openings are 
equal, r^ =: r and that the scatterers have equal transmission probability tj =: t. This 
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makes the channels comparable in the ability to transmit particles. In this case the 
merit of efficiency of our engine is simplified to 

Imagine now that the scatterers act as energy window ffiters. Then by narrowing the 
window, the variances Vi will decrease and the averages Si will converge to some definite 



value. If the numerator in (26 ) does not vanish, then the merit of efficiency y will diverge 
and consequently 

lim 7]* = ?7Carnot • (27) 

We conclude that efficiency of heat engines or refrigerators based on our model can 
be improved by strongly narrowing the energy window of transmission at the price of 
smaller power output P* or heat current Q*, depending of working regime. The idea 
of using narrow energy windows to improve thermoelectric efficiency is not new [2]. 
However it is interesting that such principle has a quite general validity. 

5. Engine with a window-shape scatterers 

Here we study a specific model of an engine (heat engine and refrigerator) with two 
(n = 2) one-dimensional {d = 1) channels with equal relative effusion probabilities 
ri = r2 = ^. The particles are effused into a channels from side z/ with momentum 
p G M distributed with probability density 

Puip) = me-^^^'Oia^p) . 

In the middle of each channel there is a point-like scatterer, which allows transmission 
only of particles with kinetic energy in a specific interval and reflects all others. The 



on-shell transmission function of the scatterer ipi ( 14 ) in ith. channel has a window shape 
in the energy axis that is given by 

UW) = e{xr + Ax. - w)e{w - Xi) , 

where Xi is the beginning of the energy window and A^j is its width. This case is 
schematically depicted in figure [2j 

The dynamics of particles in channels is conservative and symmetric w.r.t. time 
reversal. Due to simplicity of dynamics the transmission function of the channel ([I]) can 
be written explicitly and reads 

Tu^i{x,p, Ui) = Ipii-p"^ - -CTuUi) . 

In the presence of electric field, particles escape from the channels in a finite time of 
order of magnitude 0((minj |?7j|)~^/^). On the same time scale the engine relaxes into 
a non-equilibrium stationary state. 
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Figure 2. Schematic picture of the engine with the channel scatterers in behaving as 
energy windows. 



The performance of such engines can be evaluated analytically but optimization, 
in general, can only be done numerically. In the following we study engines in which 
scatterers have transmission functions of the same width A^. In particular, we set 
the electric field in the second channel to zero and optimize efficiency by fine-tuning the 
electric field in the first channel. In figures [3] and |4] we show the optimal performances of 
the heat engine and the refrigerator, respectively, in the nonlinear temperature regime 
as a function of Xi fo^' two different A^. In figure [3] we plot the optimal efficiency 
?7*(Xi,X2) and the corresponding power P*{xi,X2), which both turn out to be roughly 
symmetric functions. The symmetry is broken as we apply the electric field just in 
the first channel. The optimal efficiency, near to Carnot's, is reached on an area lying 
parallel to the symmetry line Xi = X2 and, due to asymmetry in the setup, also along 
a steeper line X2/X1 = const. A high value of rj* and with simultaneous high power 
output P* can be achieved in an area near to the Xi axis. Note that the maximal 
possible efficiency slightly increases by decreasing the width Ax of the energy window 
in the scatterers. 

The refrigerator is constructed similarly as the heat engine. The only difference is 
that we now apply an electric field to revert the direction of the heat current so that 
the heat is sucked by the engine at the cold bath Q|l > (we take T^ < Tr), and 
thereby the bias voltage performs some work (P < 0). We find the optimal performance 
of the refrigerator by numerically maximizing the efficiency rj = Q\h/{—P) with the 
constrains QIl > and P < 0. This is not a trivial procedure and can easily fail, if 
the optimal electric fields are either small or large. In figure |4] we present the relative 
optimal efficiency ?7*/?7carnot and the corresponding heat current Q* as a function of 
the energy window thresholds Xi ^^ the nonlinear temperature regime. The results are 
shown for two widths Ax of energy windows. We see that the dependence of ?7*/?7carnot 
and P* on parameters Xi and Ax is very similar to the case of optimally performing 
heat engine, as expected from the linear theory. The area, where the optimization fails 
or where efficiency is negative, is painted black and is located along the X2 axis. This 
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Figure 3. The relative optimal efficiency f7*/f7carnot (a) and the corresponding rescaled 
power aP*/[r(77cainotAx)^] (b), where for convenience we take a — 13, of the heat 
engine at temperatures Tl = 1 and Tr — 2 {T = 1) and 7 = 1 for two different width 
of energy windows Sx- Here Tycarnot = h 



show that the refrigerator can not be reahzed for arbitrary setup of scatterers. 

In the hnear temperature regime we are particularly interested in the limit of narrow 
energy windows. In this limit we can analytically compute the optimal efficiency r]* and 
the corresponding power P* and heat current Q*. Let us consider scatterers with energy 
windows of the same width Ax and at mean energies Xi ^-^id X2 i^ ^^^ fi^^^ ^^^ the 
second channel, respectively. In the discussed limit Ax — )■ we may write 



-[ne-''' + r2e-''^]aAx + 0[{Ax) 



Q* 



6/3 



[rie-^i + r2e-^2]aAa; + 0[{Ax) 



rf 



Ax 
a 



(28) 
(29) 
(30) 



1 - =^ + 0[(Ax)' 

^Carnot 

where we define rescaled window thresholds Xi := f3xi, the rescaled widths Ax := (3 Ax 
and a constant 



v^l 



X2 — Xi 



2 cosh 



X2 - Xi + log( 



Ti- 



lt is clear that by decreasing the rescaled width Ax to zero, the efficiency r] (30) 
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Figure 4. Refrigerator relative optimal efBciency 77*/77carnot (a) and corresponding 
heat current AbQ* /{T-q^^^^^^Ax) (b). Here jycarnot — 1- For details see the caption 
under figure [3j 



increases linearly up to Carnot's and the leading order of P* (28) and Q* (29) decrease 
quadratically to zero. 



Such behavior of efficiency is expected and agrees with the result (27) and with the 



energy filtering mechanism discussed in the context of quantum mechanics [2J, |5j . 



6. Heat engine with magnetic fields 

In this section, we consider a slightly different setup where time reversal invariance of 
the scattering mechanism is broken. We again take a model with two channels that 
are two dimensional d = 2 straight leads with hard walls. The ith channel is divided 
into three sections of length aj, 6j and Cj (listed from left to right), see figure [s] The 
channel's length is L = Oj + 6j + q and its width is hi. The magnetic field Bi is applied 
to the middle section in the z direction and is described by the "synchrotron orbital 
frequency" 



eBi 



m 



where e and m are the charge and the mass of the particle, respectively. In addition to 
the magnetic field, we have a homogeneous electric field Ei along the channel with 
voltage Ui = —EiLi between the ends. The dynamics in the channels is energy 
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Figure 5. An illustration of trajectories in a channel with the middle section immersed 
into a transverse magnetic field. The lower panel shows the case of a trajectory effused 
from the other side. 



conserving, but not symmetric w.r.t. time reversal. We may decompose the motion 
into bounces between channel boundaries of different sections that we may think of as 
nodes in a scattering graph. Then the scattering of particles is just a way to move 
in the scattering graph. In the considered case the scattering graph is infinite, due to 
possibility of very small loops on the boundaries. However, the total length and time 
spent in these loops is proportional to the length of the boundaries and therefore is 
finite. Consequently, the particles in the presence of the electric field are ejected from 
the channels in a finite time. 

The dynamics of particles in the lead has certain symmetry. In the absence of 
electric field each trajectory starting on the left side has a counterpart starting on the 
right side obtained by reflecting the trajectory through the center of the middle section. 
This means that at Ei = and with equal bath temperatures Tl = Tr we have tL,i = ^R,i 
and gL,i = (lR,i- The transport coefficients t^^i and q^^i depend only on the absolute value 
of the magnetic field. Indeed change in the sign of the magnetic field just mirrors the 
dynamics around the middle of the channel. 

In the following we optimize the efficiency as a function of the electric field at the 
fixed orbital frequencies fi, > 0. In numerical experiments the channels are taken of 
equal width, implying ri = r2 = ^, and with the unit average injection rate 7 = 1. 



6.1. Linear temperature- difference regime 

In order to discuss the heat engine in the linear temperature difference regime using the 



linear expansion of currents (20) we have numerically checked the validity of Eqs. (15). 
We calculate the transport coefficients t] , g| and k} (13) at the mean reciprocal 
temperature /3 = 1/T using the transmission function of the channels at zero electric 
field Ti,^i\ui=o ■ The sections without the magnetic and electric field are perfect conductors 
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and therefore transmission through the ith channel is determined by the transmission 
properties of the middle section with the magnetic field of strength fij, length bi and 
width hi. Consequently the functions Tu,i\ui=o in ( 13 ) can be replaced by the transmission 



function of the middle section TB{y,Px,Py',^i,bi,hi) and we can write the transport 
coefficients as 



^,<ir\kr: 



1 

hi 



hi 



dy / dp, / dp,P{p;P)TB{y,p;n„bi,hi){l,W{p),W{p)' 



where p = {px,Py) is the momentum and W{p) = ^{p^ +Py) is the kinetic energy of the 
particles effused into the channels. By closer analysis of dynamics we find that transport 
coefficients can be conveniently expressed as 



(f(0) „(0) , (0)^ 
K'^i iHi 1 '^i J 



2/3 1 



dWW^^e-^'^i) 



W 



kA{1,W,W^' 



(31) 



with the on-shell transmission probability ip defined as 

ijj {w, k) 



- / dy ds tb [y, Vl - s^, s, y/w, k, 1 



where we introduce the energy scale Wo^i = ^{QihiY and the ratio between the length of 
the section and its width k, = bi/hi. The energy scale is proportional to the ratio between 
the critical temperature as defined in the paper [9j and the average bath temperature 
T. In general ip{w, k) is a piecewise smooth function of both parameters as we can see 
in figure ([6]). Shorter is the section in comparison to the width, higher is the average 
transmission probability. As expected, the transmission probability at given w increases 




Figure 6. On-shell transmission V'(^j '^) through a channels with homogeneous 
magnetic field as the function of momentum |p| for various k. 

with decreasing k. If the section with the magnetic field is longer than its width, then 
ip becomes independent of k, meaning 



ipi^w, n) = ip{w, 1) 



for K > 1 . 
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Additionally, ip{w, k) as a function w is independent of k in the limit of small and large 
w. The asymptotic expansion in these two limits reads 

TT 1 1 

?/'(ti;, k) ~ — w2 tf — )■ and V"!^; ft) ~ 1 — 777^;^"'^ + 0{w^'^) w ^ 00 . 

In the absence of an electric field, if) can be obtained analytically, but it is difficult. 
Nevertheless, we work here with a numerically computed if) enhanced by analytical 



approximations. By using ip to calculate coefficients (31 ) we study the optimal efficiency 
7]* and the corresponding power P* of the heat engine given by equations (22) and (25), 
respectively, for different value of ratio /t. The results are shown in figure [7} where we 
plot r]* and P* as a function of the rescaled energy scales K^(3Wo^i. Because the channels 
are identical, 77* and P* are symmetrical on exchange of parameters. The areas of high 
optimal efficiency 77* and high power P* overlap and are located along the abscissa and 
ordinate axes in PVo,i x W^o,2 parameter space. For a given k, the maximum of the 
optimal efficiency is achieved on lines (PFo,i; 0) and (0, ^^0,2)- The position of maximal 
optimal efficiency VFo,max and its value rj^^^^ depend on k, as we see in figuresjSja andjSjb. 
The maximal efficiency rj^g^^ and the relative energy scale /3VFo,max defining its position 
are increasing as we decrease k. Numerical calculations of T/j^ax ^^^ difficult to perform 
for K near to zero. In order to obtain results in this limit, we approximate the on-shell 
transmission probability with its approximate limiting form, given by 

,,,,,,[ ^w^ : w< ^{l - VI - rcKf 

I 1 — i^KW 2 : otherwise 

within numerically determined absolute error \ip{w, K,)—ipa{'w, k)\ < 5-10^^ for k, < 10^^. 
By using this approximation we see that the maximal optimal efficiency increases with 
decreasing « and converges in the limit k — )■ to its upper bound 



Vn 



^Carnot 



= 0.0373 ± 0.001 

K=0 



at which the relative energy scale is k'^I3Wq ~ 6.55 ±0.01. Numerical results in figure [8|b 
support this finding. We see that the maximal possible efficiency in the linear regime is 
below 4% of that achieved in the Carnot's cycle. 

6.2. Non-linear temperature difference regime 

In the non-linear temperature regime, where |Tl — Tr| > T, we numerically obtained 
the transport coefficients t^^^i and q^^i ^ and study the performance of the heat engine. 
In order to obtain heat engine's optimal performance at fixed setup, we numerically find 
the optimal efficiency w.r.t. the electric potential f/j. 

As an example, we present in figure [9] the optimal efficiency 77* and corresponding 
power P* as functions of relative energy scales PWo,i for i = 1,2. In the example the 
average temperature is T = 1, the electric field is applied only in the first channel and the 
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Figure 7. The relative optimal efficiency 77*/77carnot (a) and the corresponding 
rescaled power ^*/(?7carnot^) ^ ^ function of relative energy scales /3K^Wo,i at mean 
temperature T — 1 for different values of k. 



lengths of different sections and widths of the channel are equal Oj = 6j = q = /ij = 1. In 
this setup we fix Kj = 1. The rj* and P* as the functions of /3VFo,i are slightly asymmetric, 
due to application of bias voltage to the first channel only. Their functional dependence 
is not much different than that obtained in the linear temperature regime. We find that 
the relative maximal optimal efficiency ?7max/''?Camot slightly increases when we increase 
the temperature difference between the baths. By increasing the temperature difference, 
at fixed average temperature, we also increase the power and heat current since they 
are approximately proportional to r/camot and ^camot' respectively. 
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Figure 8. Relative maximal optimal efBciency '7max/'ycainot (a) and its position f3Wo 
(b) for different values of k. 
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Figure 9. The relative optimal efficiency rj* /rjcurnot (a) and the corresponding rescaled 
power P* / {ilcainot'^) ^^ ^ function of relative energy scales (3Wo^i for temperatures 
Tl = 0.55 and Tr = 1.45 and their mean T = 1. 



7. Conclusions 



We have demonstrated the operation of a thermoelectric heat engine and refrigerator 
in terms of a classical mechanical model, which combines the deterministic classical 
scattering dynamics and stochastic heat reservoirs. The advantage of our simple model 
is that it is analytically treatable and that it allows for explicit exact results. In 
particular we have been able to disclose various situations in which the efficiency can be 
close to Carnot's and thus pose a promising alternative for technological application, in 
particular at a nanoscale, where the motion between the baths could be approximated 
by dissipationless dynamics. 
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Appendix A. Calculating the classical scattering matrix 

The channels joining thermo-chemical baths are treated as two-port scatterers for 
particles of d degrees of freedom. Let us consider the ith channel connected to the 
left and to the right bath through openings in space Fl,! and Fr^j, respectively, of an 
equal area Ai^i = Ar j ='■ Ai. The transport properties of the channels are controlled 
via homogeneous electric field with the voltage Ui between the ends. We image having 
a stationary transport in the channel with particles being injected and ejected at both 
sides at constant rate. The probability density of particle's positions and momenta that 
are injected at u side is given by 

where we introduce a momentum space for particles traveling in u direction M^ = 
Mo-^ X M.'^^^. Similarly, the probability density of particles ejected at side u is 



PoutA^^P) ix,p) eV^^iX M- 



<u 
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We assume pin,u{x,p) = for ai^pi < and Pont,u{x,p) = for a^pi > 0. Using the 
introduced notation we may write the probabihty density of ingoing particles pin on the 
phase space V^n,i = (rL,j x Ml) U (Fr^j x Mr) and the probabihty density of outgoing 
particles pout on the phase space Vout,i = (Fl.j x Mji) U (Fr^j x Ml) are defined as 

Pin = Pin,L + Pin,R and Pout = Pout.L + Pout.R , 

respectively. The two-port classical scattering operator S associated to the channel is a 
linear map of pin into pout^ 

Pout = '-'Pin • 

The operator S depend on dynamical properties of the channel. By using projectors P„ 
defined as P^f = 0{(7^pi)f, we can decompose the operator S into four operators: 

Th = PhSPh , -Rl = PrSPl , ^R = -Pr'S'Pr , -Rr = Pl'S'Pr . 

These operators describe (listed in the same order as above from left to right) the 
transmission from the left to the right side, the reflection from the left back to the 
left side, the transmission from the right to the left side and the reflection from the 
right back to the right side. Notice the sum of these operators is again the operator 
S. The dynamics is deterministic and pointwise. Consequently, the trajectory's entry 
point z = {x,p) can be uniquely connected to the exit point z' = {x',p') via a map 
z' = (f){z) and S" is a Koopman operator. This means S6'^ = (^^(z'), where we use Dirac 
delta S'^{z) := 6{z — z'). By using these facts we express the transmission function Tj^^j 
Q as 



T,,{x',p', U,) = / dz {TJ,,){z) = / dz {S6,>){z) . (A.l) 

The transmission function Tj^i of a channel is an essential ingredient of the presented 



theory of heat engines and refrigerators, and relation (A.l) connects it with a scattering 
operator S. As we will show in the following, the operator S can be systematically 
obtained for structurally complicated channels, when scattering in its parts is known. 

Let us assume that a channel can be divided into two parts with a common cross- 
section and simpler dynamics, which can be described by the scattering operators 5*1 and 
5*2. The operator Si, can be decomposed into transmission and reflection operators Tj ,^ 
and Ri^u- By knowing them, we can construct the transmission and reflection operators 
Tiy and R^ corresponding to the scattering operator S of the whole channel via the 
concatenation formulas 

Rl = Rl,L + TiYiR2,LL Ti L , 7l = T2_L-Z^ ^1 L , 
-Rr = -R2,R + ^2,L-Rl,R-^' ^2,R , ^R = Ti^L' T2^r , 

with L = 1 — -Ri,R-R2,L and L' = 1 — R2,lRi,r- This is analogous to concatenation of 
quantum scattering matrices discussed in [10]. Note that the inverse operator (1 — A)^^ 
is just an symbolic abbreviation for the series ^^q^*- ^^ order to simplify the writing 
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we represent the concatenation of two scattering operators 5*1 and 5*2 into a resulting 
scattering operator S* by a non-linear and non-commutative product labeled by the 
symbol © yielding 

O ^ Ol 02 • 

For example, let us consider a straight empty one dimensional channel of length L in 
the middle of which is a point scatterer that maps the incoming particle momentum p 
into p' = 0{e)p, where O orthogonal matrix depending on the energy e = ^p'^p. Such 
maps obeys the energy conservation. The scattering operator associated with the empty 
wires of length |L with the voltage ^U between ends is defined as 

{E(f)){pi,p2, . . . ,Pd) = 0(-a)0(-pi,p2, • • • ,Pd) + ^ ^ 0(sign(pi)Va,p2, ■ ■ • ,Pd) 

with a = p\^ sign(pi)f/. The scattering operator corresponding to the scatterer in the 
middle of the channel is written as 

(P0)(p) = 0(O(e)-V). 

Then we can write the scattering operator corresponding to the whole channel as a 
product 

S = EQ)P(dE. 

The latter can be calculated systematically. Such kind of channels are used for 
conduction of one-dimensional particles, where 0(e) G { — 1,1}, in a realization of a 
heat engine discussed in section [5] 

Appendix B. Relaxation times 

We consider a heat engine or refrigerator with multiple channels and no particles initially 
in the channels. Then at time t = 0, the particles start to be injected into the channels 
from the left and right bath with the injection rates 7l and 7pj,, respectively. The 
number of particles in channels A^(t) increases with time. Eventually it saturates to 
a stationary value and so the heat engine enters into a non-equilibrium stationary 
state. The characteristic time scales of this relaxation process are essential for practical 
implementations of heat engines and refrigerators. These can be obtained through the 
study of functions iV(t), defined more precisely in the following. 

Let us introduce time function 0,^^j, where (pi,,i{q,p) > represents the time needed 
for a particle entering ith channel from side z/ at the point (g, p) G Vu^i = T^^i x M^ to 
return into the baths. The time function can be measured for channels and in some 
simple cases even analytically expressed. The number of particles in ith channel coming 
form side u is given by 

N^,i{t) = -— dq dpP{p,au,^u)TuAl^P^Ui)mm{4)^^^{q,p),t) . (B.l) 
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Then the total number of particles in the channels is equal to 

Nit) = J2niuN,,,{t). 

In order to obtain the relaxation time scales we need to understand asymptotics of 
N{t), which is determined by effective injection rates rj7j, and asymptotics of individual 
contributions Nj^^i. It is difficult to say something definite about the latter. But if 0j,^j 
is continuous except on Uu^i C Vi,^i of zero-measure then Ni,^i is also continuous. From 



definition (B.l) it is clear, that N{t) is non-constant in the limit t — )■ oo, if 0j,j has a 
singularity on subset oiUu,i- This nontrivial asymptotics of N{t) is determined by the 
functional dependence of (p^^i in the vicinity oiU^i. 

As an illustration, let us consider an empty one-dimensional channel of length L 
with one- dimensional particles injected from one side at the rate 7. The time function 



of such channel depends on the voltage U and reads 

2L 



u < 



0{V^-\p\Y-^ + 0{\p\-^)^== : otherwise ' 

v+yp +u 

with u = 2 sign(p)f/. By using we can compute the number of particles in the channel 
coming form one side Ni{t) as a function of time t. In the presence of the electric field 



is bounded from above by to = L/ ^J2\U\, which represents the time needed to achieve 
the stationary state. The stationary number of particles is A^o '■= ^™^t-¥oo ^lif) oc 7^0- In 
the absence of the field (f){p) does not have an upper bound and behaves as 0(p) x p~^ 
near p = 0. Consequently, the number of particles Ni converges slowly in algebraic 
manner as 



N,{t) - No + ^ hh- —^ ] +0{t' 



with the stationary number A'"o = '~^L\f^Txj2 and the time scale of relaxation t\ = 
La//3/2. In view of these results, we expect in general for channels composed of parallel 
walls that in the absence of the field the convergence towards the stationary state is 
algebraic. Therefore to achieve the stationary state in a finite time in such cases it is 
practical to apply electric fields in all channels. 



